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ABSTRACT 

PSR J1119-6127 is a high magnetic field (B=4.1xl0 13 Gauss), young (<1,700 
year old), and slow (P=408 ms) radio pulsar associated with the supernova rem- 
nant (SNR) G292.2-0.5. In 2003, Chandra allowed the detection of the X-ray 
counterpart of the radio pulsar, and provided the first evidence for a compact 
and faint pulsar wind nebula (PWN). We here present new Chandra observations 
which allowed for the first time an imaging and spectroscopic study of the pulsar 
and PWN independently of each other. The PWN is only evident in the hard 
band (above ~2 keV) and consists of jet-like structures extending to at least 7" 
from the pulsar, with the southern 'jet' being longer than the northern 'jet'. The 
spectrum of the PWN is described by a power law with a photon index r~l.l for 
the compact PWN and ~1.4 for the southern long jet (at a fixed column density 
N H =1.8xl0 22 cm" 2 ), and a total luminosity 2^(0.5-7.0 keV)~4xl0 32 ergs s -1 , 
at an assumed distance of 8.4 kpc. The pulsar's spectrum is clearly softer than 
the PWN's spectrum. We rule out a single blackbody model for the pulsar, and 
present the first evidence of non-thermal (presumably magnet ospheric) emission 
that dominates above ~3 keV. A two-component model consisting of a power law 
component (with photon index 5-2.0) plus a thermal component provides 
the best fit. The thermal component accounts for the soft emission and can be fit 
by either a blackbody model with a temperature fcT~0.21 keV, or a neutron star 
atmospheric model with a temperature &T~0.14 keV. The efficiency of the pulsar 
in converting its rotational power, E, into non-thermal X-ray emission from the 
pulsar and PWN is ?»5xl0~ 4 , comparable to other rotation-powered pulsars with 
a similar E. We discuss our results in the context of the X-ray manifestation of 
high-magnetic field radio pulsars in comparison with rotation-powered pulsars 
and magnet ars. 

Subject headings: ISM: individual (G292.2— 0.5) - - pulsars: individual (PSR 
J1119— 6127) — supernova remnants — X-rays: ISM 
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1. Introduction 

A pulsar wind nebula (PWN) is a bubble of relativistic particles and magnetic field, 
powered by the rotational energy of a young energetic pulsar. The PWN's non-thermal 
emission arises from synchrotron radiation of the high-energy particles injected by the pulsar 
and accelerated at the site where the relativistic pulsar wind is decelerated by the confining 
medium. As such, PWNs represent an ideal laboratory to study the physics of neutron stars, 
particle acceleration, and relativistic shocks. 

The Crab nebula has been viewed for almost three decades as the canonical example 
of how a fast rotation-powered pulsar dumps its energy into its surroundings. However, 
the past several years have witnessed a synergy of high-energy observations pointing to a 
growing class of compact objects which manifest themselves differently from the Crab. These 
include the anomalous X-ray pulsars (AXPs) and the soft gamma-ray repeaters (SGRs), 
commonly believed to be highly magnetized neutron stars, or magnetars (S~10 14 -10 15 G), 
powered by magnetic field decay rather than by rotation, and characterized by relatively 
slow periods (5-12 seconds; see Woods and Thompson 2006 for a review). So far, there is no 
evidence of PWNs surrounding these objects. Furthermore, they were believed to be radio 
quiet, however recently radio emission has been observed from the transient magnetar XTE 
J1810— 197 (Halpern et al. 2005, Camilo et al. 2006) and 2-second radio pulsations have 
been discovered from a new magnetar (IE 1547.0—5408, Camilo et al. 2007). 

Interestingly, there are about half a dozen radio pulsars with spin and magnetic field 
properties intermediate between the 'classical' Crab-like pulsars! and the magnetars: their 
spin periods range from a fraction of a second to ~7 seconds, and their magnetic fielc0 
is close to the QED value of 4.4 xlO 13 Gauss (the field at which the electron's rest mass 
energy equals its cyclotron energy). The question whether these pulsars should display more 
magnetar-like or Crab-like properties is still under debate (see e.g. the references in Table 3). 
X-ray observations targeted to study these objects offer the best window to address this 
question and better understand their emission properties. The search for and study of any 
associated PWN should in particular shed light on the nature of their high-energy emission 
in comparison with the other pulsars. 

PSR J1119— 6127 is a member of this small, yet growing, class of objects. It was 
discovered in the Parkes multi-beam pulsar survey (Camilo et al. 2000). It has a rotation 



3 A Crab-like pulsar here refers to a rotation-powered pulsar powering a PWN 

4 the magnetic field is derived assuming magnetic dipole radiation and is related to the spin period, P, 
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period P of 408 ms, a large period derivative P of 4x10 12 , a characteristic age ^ of 
1600 yrs, a surface dipole magnetic field B of 4.1 xlO 13 G, a spin-down luminosity of 
2.3 xlO 36 ergs s -1 , and a braking index of 2.9 implying an upper limit on its age of 1,700 years 
(under standard assumptions of simple spindown). No radio PWN was detected around the 
pulsar despite its youth (Crawford et al. 2001). The pulsar was found within the supernova 
remnant (SNR) G292.2— 0.5 of 15' in diameter, discovered with the Australia Telescope 
Compact Array (ATCA) (Crawford et al. 2001), and subsequently detected with X-ray 
observations acquired with the ROSAT satellite and the Advanced Satellite for Cosmology 
and Astrophysics, ASCA (Pivovaroff et al. 2001). 

The X-ray counterpart to the radio pulsar was first resolved with Chandra, which also 
revealed for the first time evidence for a compact (3"x6") PWN (Gonzalez & Safi-Harb, 
2003; hereafter GSH03). The pulsar was subsequently detected with XMM-Newton and was 
found to exhibit thermal emission with an unusually high pulsed fraction of 74±14% in the 
0.5-2.0 keV range (Gonzalez et al. 2005). The XMM- Newton-derived temperature was also 
unusually high, in comparison with other pulsars of similar age or even younger, like the 
pulsar in the SNR 3C 58 (Slane et al. 2002). 

While a spectrum of the pulsar has been obtained with the previous Chandra and 
XMM-Newton observations, the Chandra spectrum suffered from poor statistics and it 
wasn't possible to extract a spectrum from the PWN; and the XMM-Newton spectrum of 
the pulsar was contaminated by the PWN's spectrum due to the large point spread function 
(PSF) of XMM-Newton. 

In this paper, we present new observations of the pulsar obtained with Chandra and 
targeted to constrain the properties of the pulsar and associated PWN. In a follow-up paper, 
we present the study of the associated SNR, to complement the previous study of the western 
side of the remnant by Gonzalez & Safi-Harb (2005). 

The distance to the pulsar was determined to be in the 2.4-8.0 kpc range (Camilo et al. 
2000). Using the fitted X-ray column density to the pulsar and that to its associated SNR, 
GSH03 estimated a distance in the 4.0-12.6 kpc range, but the location of the source with 
respect to the Carina spiral arm (Camilo et al. 2000) puts an upper limit on the distance of 
~8 kpc. The Cordes-Lazio Galactic free electron density model (NE2001, Cordes & Lazio 
2002) suggests a much larger distance of 16.7 kpc, while Caswell et al. (2004) estimate a 
distance of 8.4 kpc using neutral hydrogen absorption measurements to the SNR. In this 
work, we adopt the distance to the SNR of 8.4 kpc, and scale all derived quantities in units 



The paper is organized as follows. In §2, we summarize the Chandra observations. In 
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§§3 and 4, we present our imaging and spectroscopic results, respectively. Finally, in §5 we 
present our discussion and conclusions. 

2. Observations and Data Analysis 

PSR J1119— 6127 was observed with Chandra in cycle 5 on 2004 November 2-3 (obsid 
6153), and October 31-November 1 (obsid 4676). We combine these new data with the ob- 
servation conducted during cycle 3 (obsid 2833, GSH03) in order to improve on the statistics. 
All three observations were positioned at the aimpoint of the back-illuminated S3 chip of the 
Advanced CCD Imaging Spectrometer (ACIS). X-ray photon events were acquired in the 
timed-exposure readouts, at a CCD temperature of — 120°C and a frame readout time of 3.2 
sec, and telemetered to the ground in the 'Very Faint' mode. The roll angles were chosen 
such as the AO-3 data covered the western side of the SNR and the AO-5 data covered the 
eastern side (see Fig. 1). The data were then reduced using standard CIAO 4.0 routine^. We 
have created a cleaned events file by retaining events with the standard ASCA grades 02346 
and rejecting hot pixels. As well, the data were filtered for good time intervals, removing 
periods with high background rates or unstable aspect. The resulting effective exposure time 
was 56.8 ksec for obsid 2833, 60.5 ksec for obsid 4676, and 18.9 ksec for obsid 6153, yielding 
a total effective exposure time for the 3 combined observations covering the field around the 
pulsar of 136.2 ksec. 

3. Imaging results 
3.1. G292.2-0.5 

Figure 1 shows the Chandra cycle 3 (covering the western side of the remnant) and 
cycle 5 (covering the eastern side of the remnant) observations overlayed on the ASCA 
hard band (2-10 keV) image, and with the radio extent of the SNR shown as a black circle 
overlayed. The central circle denotes the position of the pulsar. The image illustrates 
the power of Chandra in resolving the X-ray emission and shows the overlap between the 
Chandra observations leading to a deep exposure of the pulsar and central regions of the 
SNR. The pulsar and surrounding region will be the focus of this paper. The properties of 
the diffuse emission from the SNR will be presented in a separate paper. 



5 http:/ /cxc. harvard.edu/ciao 
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Figure 2 shows the combined Chandra tri-color image, zoomed on the pulsar region, 
and obtained as follows: the data were firstly divided into individual images in the soft (0.5- 
1.15 keV, red), medium (1.15-2.3 keV, green), and hard (2.3-10.0 keV, blue) bands (the 
energy boundaries were chosen to match those in GSH03). Individual background images 
were then obtained from the blank-sky data sets available in CIAO. For each observation, 
background images were reprojected to match the corresponding events file. The background 
images were then divided into the same energy bands mentioned above, and then subtracted 
from their corresponding source+background image. The resulting images were subsequently 
adaptively smoothed^] using a Gaussian with a—1" for significance of detection >5 and up 
to a=5" for significance down to 2. The individual background-subtracted images were 
then combined to produce the image shown in Figure 2. As in our previous X-ray study 
(GSH03), the X-ray emission can be resolved into several X-ray sources (tabulated in Table 2 
of Gonzalez & Safi-Harb 2005), surrounded by diffuse emission from the SNR interior. As 
well, the pulsar and associated PWN immediately stand out as a hard (blue in color) source. 



3.2. PSR J1119-6127 and its PWN 

The peak of the X-ray emission from the pulsar is positioned at: aj2ooo = H /l 19 m 14 f 526, 
5 J2000 = -61°27'49"3, which match within error the radio coordinates of the pulsar: ctj2ooo 
= ll h 19 m U 8 .3 and 5j 2000 = -61°27'49"5 (with 0.3" error). 

In order to compare the source's spatial characteristics with Chandra's PSF, we per- 
formed a two-dimensional spatial fit to the new Chandra data using the SHERPA software^. 
We created images of the source in the soft (0.5-1.15 keV), medium (1.15-2.3 keV), and hard 
(2.3-10 keV) bands. We subsequently created normalized PSF images at the energies charac- 
teristic of the source's energy histogram. These images were then used as convolution kernels 
when fitting. The soft band image yielded a FWHM value of 0.8" consistent with a point 
source, the medium band image gave a slightly higher value of 0.9", while the hard band 
was best described by an elliptical Gaussian function with FWHM 0.9"xl.2", confirming the 
extended nature of the PWN in the hard band, a result that is consistent with GSH03. 

In Figures 3 and 4, we show zoomed in images of the PWN. Figure 3 shows the raw soft 
band (0.5-2 keV) and hard band (2-10 keV) images centered on the pulsar, and Figure 4 



6 the 'csmooth' CIAO tool is based on the 'asmooth' algorithm (Ebeling et al. 2006) which has the 
advantage of visually enhancing extended features in noise-dominated regions. We note that the compact 
PWN and the southern jet display a similar morphology using a simple Gaussian smoothing. 

7 http:/ /cxc. harvard.edu/sherpa/ 
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shows the same images adaptively smoothed using a Gaussian with a—1" for a significance 
of detection greater than 5 and up to a=3" for a significance of detection down to 2. As 
mentioned earlier, the PWN structures (further discussed below) are also visible using a 
simple Gaussian smoothing. We verified that they are also evident in the hard band images 
for the individual observations. Furthermore, using the celldetect CIAO routine, we did 
not find any point sources within the PWN, thus ruling out the contamination of the PWN 
structures by background sources. 

The images clearly show that the PWN is harder than the pulsar and can be resolved 
into elongated 'jet'-like structures extending to at least 7" north and south of the pulsar. 
This confirms the evidence of a hard PWN associated with this high-magnetic field pulsar 
(GSH03) but also shows that the PWN is larger than originally thought. The PWN is 
elongated in the north-south direction, with the northern extension (which we refer to as 
a 'jet') appearing tilted to ~30° west of north, and the southern jet being more elongated 
and ~2x brighter than the northern 'jet'. In Figure 5, we highlight the first evidence 
for a more extended (>20") faint and hard jet-like structure southwest of the pulsar. We 
accumulate 100±23 background-subtracted source counts from this region, which on top of 
99±21 background counts (normalized to the source area), imply a ~7±2cr detection (0.5-7.0 
keV). It is not clear whether this structure is an extension of the southern compact jet seen 
closer to the pulsar or an SNR filament. However, its morphology (in particular its knotty 
structure which is similar to Vela's long jet, Kargaltsev et al. 2003) and spectrum (see §4.2) 
suggest an association with PSR J1119— 6127's PWN, a result that can be confirmed with 
a deeper exposure. We note that while the XMM-Newton MOS2 image shows a hint of 
a southern extension, the large and distorted PSF of the MOS2 camera does not allow a 
confirmation of this jet (M. Gonzalez 2008, private communication). 

In the following, we refer to the PWN region shown in Figure 4, right (size ~6" in the 
east-west direction by ~15" in the north-south direction ) as the 'compact PWN', and to 
the extended jet-like structure southwest of the pulsar (size ~6"x20", see Fig. 5) as the 
'southern jet'. We note here that it's possible that the compact PWN corresponds to a torus 
viewed edge-on. However, the southern extension (if proven to be a jet associated with the 
pulsar) would falsify this interpretation as the jets are expected to be perpendicular to the 
torus. We further discuss these features in §4.2 and §5.2. 
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4. Spatially resolved spectroscopy 

The spectral analysis was performed using the X-ray spectral fitting package, XSPEC 
vl2.4.C@, and constrained to the 0.5-7.0 keV band as the spectra provided little useful infor- 
mation outside this range. The pulsar's spectrum was grouped with a minimum of 20 counts 
per bin, while the PWN's spectrum was binned using a minimum of 10 counts per bin. 
Errors are at the 90% confidence level throughout the paper. 

4.1. PSR J1119-6127 

The spectrum of PSR J1119— 6127 was extracted using a circle centered at the posi- 
tion of the pulsar (§3) and with a radius of 2.5", thus encompassing 90% of the encircled 
energj^]. The background was extracted from a ring centered at the pulsar, and extending 
from 3" to 4". This has the advantage of minimizing contamination by the PWN. We note 
however that the PWN contamination is negligible as its surface brightness is nearly two 
orders of magnitude smaller than the pulsar's, and so the spectral results are practically 
independent of the background subtraction region (see e.g. Safi-Harb 2008 for a preliminary 
analysis of the pulsar's spectrum using a source-free background region taken from the S3 
chip 10 !) . The background-subtracted pulsar count rate is (5.2±0.3) xl0~ 3 counts s _1 for the 
2833 observation, (5.8±0.3) x 10~ 3 counts s" 1 for the 4676 observation, and (3.7±0.5) xl0~ 3 
counts s _1 for the 6153 observation, thus accumulating ~715±27 total counts from the pul- 
sar. To model the emission from the pulsar, we used a blackbody model to account for any 
surface thermal emission from the neutron star and a power law to account for any magne- 
tospheric emission. A power law model yields a column density Nh={1-2-q',3) xlO 22 cm -2 , 
a steep photon index r=3.0±0.5, and an unabsorbed flux of (1.7^q; 3 ) x 10~ 13 ergs cm -2 s _1 , 
with a reduced chi-squared value x 2 = 1.356 (i>=31 dof). A blackbody model does not yield 
an acceptable fit (x 2 =2.166, u=31 dof); and as shown in Figure 6, it does not account 
for the emission above ~3 keV. Adding a power law component to the blackbody model 
improves the fit significantly, and yields an excellent fit with the following parameters: 
AT^=(i.8+^)xl0 22 cm- 2 ,r=1.9+^, £;T=0.21±0.01 keV, xl=0M8 (i/=29 dof), with the lu- 
minosity of the power law component representing ~30% of the luminosity of the blackbody 

8 http://xspec. gsfc.nasa.gov 

9 http:/ /cxc. harvard.edu/proposer/POG/html/ ACIS.html 

10 The tabulated values in Safi-Harb (2008) correspond to binning the data with a minimum of 10 counts 
per bin, whereas the data here are binned by a minimum of 20 counts per bin. 
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component (or ~23% of the total pulsar's luminosity) in the 0.5-7.0 keV range. Figure 7 
shows the blackbody+power law model fit, and Table 1 summarizes the best fit parameters. 

The thermal and soft component is equally well fit with a neutron star magnetized 
(B=10 13 G) atmospheric model (the NSA code in XSPEC, Zavlin et al. 1996). The pa- 
rameters for the blackbody and atmospheric model are however different, and as expected 
the NSA model yields a lower temperature than the blackbody model. Both models require 
a hard photon index (r~l. 5-2.0) for the power law component. As shown in Table 1, the 
non-thermal unabsorbed flux amounts to ~30-35% of the thermal flux. In the NSA model, 
the distance was fixed to 8.4 kpc. We attempted to fit for this distance, however we find 
that it is poorly constrained (ranging from 2.6 to 30 kpc). The best fit parameters are: 
iV H =1.6xl0 22 cm" 2 , r PL =1.5, T NSA =l.Qxl0 6 K, and D=8A kpc. Therefore we restrict our 
discussion below to the fit with the fixed and reasonable distance of 8.4 kpc. 



4.2. The PWN 

The compact PWN's spectrum was extracted from the elongated nebula (~6"xl5" )(see 
Fig. 4, right) with the pulsar region removed. The background was extracted from a nearby 
source-free region. The total number of background-subtracted counts accumulated from 
the PWN is only ~103±17 counts in the 0.5-7.0 keV range, therefore restricting our ability 
to constrain the spectral parameters or to perform a spatially resolved spectroscopy of the 
northern and southern extensions. However, fixing the column density to 1.8 xlO 22 cm" 2 
(the best fit value for the pulsar in the BB+PL fit, Table 1), we find that a power law model 
yields a hard photon index r=l.llQy, x 2 = 1.23 (z/=17 dof), and a luminosity L^O. 5-7.0 
keV)~(1.6+^)xl0 32 d\ A ergs s -1 . Allowing N H to vary between (1.2-3.3) xlO 22 cm" 2 (the 
range allowed by the PL+BB model fit to the pulsar), the PWN's index is r=0.6-2.9, typical 
of PWNs. 

We subsequently extracted a spectrum from the southern 'jet' (region shown in Fig. 5) 
and used the same background region as for the compact PWN. As mentioned earlier, we 
accumulated a total of ~100±23 source counts from this region in the 0.5-7.0 keV range, 
again too small to constrain the parameters or perform a spatially resolved spectroscopy 
along this feature. Fixing the column density to 1.8xl0 22 cm" 2 , the power law model 
yields a photon index of r=1.4+g;g, x£=1.16 [y=\1 dof), and a luminosity of 1^(0.5-7.0 
keV)~(2.1± 3 ; 3 )xl0 33 d\ A ergs s -1 . These values are not unusual for PWN structures. A 
thermal bremsstrahlung model is rejected based on the unrealistically high temperature 
(A;T>10 keV). Experimenting with two-component models was not possible given the insuf- 
ficient number of counts. In Table 2, we summarize the best fit power law model results. 
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5. Discussion 
5.1. The pulsar 

For the first time, we have resolved the pulsar's spectrum and determined its properties 
independently of the PWN. We note here that the hard X-ray component in the XMM- 
Newton spectrum has been attributed entirely to the emission from the PWN (Gonzalez et 
al. 2005). We believe that this component is contaminated by the hard X-ray emission from 
the pulsar as XMM's PSF does not allow resolving both components. In fact, an on-off pulse 
spectroscopic analysis of the JMM's spectrum shows evidence of pulsed power law emission 
at the 2cr level (M. Gonzalez 2008, private communication). 

The single blackbody model is rejected (see §4.1 and Fig. 6). As shown in Table 1, the 
best fit model requires two components: a thermal component described by a BB or NSA 
model, and a power law component with a photon index 5-2.0. The column density 
is model dependent. GSH03 derived an ^=(9+3) xlO 21 cm -2 for the pulsar. Using XMM, 
Gonzalez et al. (2005) derived a higher Nh in the (1.3-2.7) xlO 22 cm -2 range (depending on 
the model used). Our best fit two-component models shown in Table 1, which adequately 
account for the hard (above 3 keV) emission, yield a column density of (1.4-1.8) xlO 22 cm -2 
in the NSA+PL model, and (1.2-3.3) xlO 22 cm -2 in the PL+BB model. These ranges 
bracket the value derived using; HEASARC's Nh too0 and overlap with the values derived 
previously with Chandra and XMM. 

Using the thermal model in the two-component fits tabulated in Table 1, we now esti- 
mate the thermal properties of the pulsar. In the blackbody fit, the radius is inferred from 
equating the isotropic luminosity (4nD 2 F, where D is the distance to the pulsar and F is 
the unabsorbed flux, tabulated in Table 2) with the surface emission from the neutron star 
(47ri? 2 o"T e jj 4 , where a is the Stefan Boltzmann constant), yielding i?=2.7±0.7 km, suggesting 
polar cap emission with a temperature of T e jj=0.21±0.01 keV=(2.46±0.12)xl0 6 K. The in- 
ferred radius is however larger than the conventional polar cap radius: ^(^-0.5 (^) 1/2 =0.25 
km (for P=408 ms). The pulsar's thermal luminosity (assuming isotropic emission) is L x (0.5- 
7.0 keV) = (1.87^)xl0 33 4 4 ergs s" 1 . 

In the atmospheric (NSA) model, the temperature is T e jj=(1.59^Q 23) xlO 6 K, implying 
a temperature as seen by an observer of T^j=g r T e ff =(1.22^^1) xlO 6 K; where g r is the 
redshift parameter given by g r =(l-2.952x ^) 0,5 =0.766 for a neutron star of mass M=1.4 



11 this tool calculates the weighted average measurement of the total Galactic column density towards the 
source's coordinates based on the Dickey & Lockman (1990) H I map. 



-10- 



solar masses and radius R=10 km. The pulsar's bolometric luminosity as seen by an observer 
and assuming thermal emission from the entire neutron star's surface is 00 (0.5-7.0 keV) 
= g 2 x(47ri?VT e 4 // ) = (2.67±1.9)xl0 33 (i 2 i4 ergs s _1 . These estimates differ from the values 
inferred by Gonzalez et al. (2005), mainly because the distance was fixed here at 8.4 kpc, 
while it was fit for with XMM. Our values for both the blackbody and NSA model are 
however consistent with those derived independently by Zavlin (2007a). 

The thermal emission from the pulsar was already discussed by Gonzalez et al. (2005) 
and Zavlin (2007a) who point out that the large pulsed fraction detected in the soft band 
is hard to interpret in either the blackbody or the NSA model (assuming uniform surface 
emission). The large pulsed fraction indicates that the thermal emission is intrinsically 
anisotropic and that the effective temperature derived in the NSA model should be used 
as the mean surface temperature (Zavlin et al. 2007b). We here point out to a puzzle 
related to the interpretation of the hard non-thermal pulsar emission, seen for the first time 
with the Chandra observations presented here. The non-thermal emission is presumably 
magnetospheric emission as seen in other young active pulsars. Such an emission is expected 
to be pulsed. As mentioned earlier, XMM could not resolve the pulsar's emission from the 
PWN, however an on-off pulse spectroscopic study of the XMM spectrum shows a hint of 
pulsed power law emission from the pulsar. A long observation with Chandra in the timing 
mode and/or a deeper XMM-Newton observation will confirm and constrain the pulsed 
component in the hard band. 



5.2. The PWN 

The PWN is well described by a power law model with a hard photon index ~1.1 for 
the compact PWN, and ~1.4 for the southern extended 'jet' (for a column density fixed at 
1.8 xlO 22 cm -2 ), and with a total X-ray luminosity of ^4xl0 32 D| 4 ergs s" 1 (0.5-7 keV). 
This represents <20% of the pulsar's luminosity, and <0.02% of E, implying an efficiency 
which is comparable to PSR B1800— 21 and other Vela-like pulsars (Kargaltsev et al. 2007). 

PWNs have been observed to have axisymmetric morphologies, including a toroidal 
structure and jets along the pulsar's spin axis. The emission from the torus is commonly 
attributed to synchrotron emission from the relativistic pulsar wind in the presence of a 
toroidal magnetic field in the equatorial plane and downstream of the shock. Jet-like struc- 
tures on the other hand are attributed to magnetically collimated winds along the pulsar's 
rotation axis. The collimation and acceleration of the jets is strongly dependent on the 
pulsar's wind magnetization a, defined as the ratio of electromagnetic to particles flux. At 
higher magnetization, equipartition is reached in the close vicinity of the termination shock, 
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and most of the plasma is diverted and collimated into a jet along the polar axis. Therefore 
the morphology of the PWN sheds light on the geometry of the pulsar and the intrinsic 
parameters of the pulsar wind (see e.g. Gaensler & Slane 2006 and Bucciantini 2008 for 
reviews) . 

The jet-like structures seen north and south of PSR J1119— 6127 could be interpreted 
as collimated outflows along the pulsar's polar axis, in which case a torus could be expected 
in the east-west direction. A torus-like structure is however barely resolved and the PWN 
emission is mostly evident in the form of elongated jets (see Figs. 3-5), suggesting a high- 
er wind. In the following, we use the spectral properties of the PWN to infer its intrinsic 
properties such as the termination shock, magnetic field, and a. 

The termination shock site (r s ) is estimated by equating the thrust of the pulsar, 47n ^ en , 

B 2 a 

with the pressure in the nebula, Prs j 2 -; where f2 is the pulsar's wind solid angle which ac- 
counts for a non-isotropic wind, and B n is the nebular magnetic field which can be approxi- 
mated from its equipartition value B eq . To estimate B eq , we use the non-thermal properties 
of the compact PWN (see Table 3), and assume the emission region to be a cylinder with 
radius r n <3"=0.12 d 8A pc and length of /~15"=0.61 d 8A pc, resulting in an emitting vol- 
ume of U~8.4xl0 53 fd\ A cm 3 ; where / is the volume filling factor. Using a luminosity of 
Lj; (0.5-7.0 keV)~2xl0 32 d\ A ergs s -1 and assuming a ratio of 100 for the baryon to electron 
energy density, we estimate a magnetic field of B eq >A.5x 10~ 5 (/d 8 . 4 )~ 2 / 7 Gauss. The corre- 
sponding nebular pressure is 1.6xl0 -10 {fd 8 A)~ 4 ^ 7 ergs cm -3 and the termination shock is 
at a radius r s ~0.06 (/d 8 . 4 ) 2/7 ^~ 1/2 pc (or 1.5" at 8.4 kpc), comparable to B1800-21 other 
PWNs (Kargaltsev et al. 2007). 

We can further check the above estimate by using the thermal pressure from the inte- 
rior of the remnant (^2n e kT) as the pressure confining the PWN. Using the kT and n e t 
values derived from our spectral fits to the SNR interior (to be presented in our follow up 
SNR paper), we derive an SNR internal thermal pressure of ~2.9xl0 -9 ergs cm -3 . The 
resulting shock radius is r s ~0.015 pc, smaller than our estimate above, but consistent for 
(/d 8 .4) 2 / 7 fi -1 / 2 =0.25. At a distance of 8.4 kpc, a termination shock radius of 0.015-0.06 
pc translates to a radius of r s =0 / .'4-l'. / 5 d 8 .4, which explains the lack of resolvable wisp-like 
structures as seen in other (nearby) PWNs. 

Using the basic Kennel & Coroniti (1984) model developed for the Crab nebula, we can 
further infer a rough value for the magnetization parameter ~0. 02-0. 25. This value 

is larger than that inferred for Crab-like PWNs (see e.g. Petre et al. 2007), but supports 
the picture for strong collimation of jets. We caution however that this estimate should not 
be taken at face value and should be only considered as evidence for a highly magnetized 
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wind, as it was inferred using the simple model of Kennel & Coroniti (1984) which assumes 
an isotropic energy flux in the wind. More recent calculations however show that if hoop 
stresses are at work in the mildly relativistic flow, jet collimation can occur in the post shock 
region and estimates of a using the Kennel and Coroniti model can be overestimated (see 
e.g. Bucciantini 2008 and references therein for a review). 

The corresponding synchrotron lifetime of a photon with an energy of E(keV) is roughly 
Tsyn^SGB 3 /^ E 1 / 2 dg^ 7 yr, where -B„4 is the nebular magnetic field in units of 10~ 4 Gauss; 
yielding an average velocity for the high-energy electrons to reach the edge of the southern 
jet of ~0.093c B_ Z J 2 E~ 1 / 2 d^J. This velocity is not unusual and is smaller than the typical 
shocked wind velocity of ~0.3c. 

We interpreted the elongated structures to the north and south of the pulsar as jets 
that could be aligned with the rotation axis of the pulsar. The northern extension is shorter 
in length and is approximately twice as faint as the southern jet. Assuming that these 
elongated features represent an approaching (south) and receding (north) jet with similar 
intrinsic surface brightness and outflow velocities, we can use relativistic Doppler boosting 
to account for the observed difference in brightness (Mirabel & Rodriguez 1999). For the 
observed brightness ratio of ~2 and an energy spectral index — a~0. 1-0.4 (Table 2), assuming 
a continuous jet, we estimate a velocity v cos #~0. 14-0. 16c, where 9 is the inclination angle 
of the flow axis with respect to our line of sight. Since cos#<l for all inclination angles, this 
represents a lower limit for the intrinsic velocity of these features, a result that is consistent 
with our estimate of the average velocity above. 

Elongated jets have been observed in other PWNs, e.g. PSR B1509— 58 powers a 4'-long 
jet (6 pc at a distance of 5.2 kpc), and Vela powers a 100"-long jet (0.14 pc at a distance of 
300 pc). Furthermore, velocities of >0.2c have been also inferred from X-ray studies of these 
jets: Vj>0.2c for B1509-58 (Gaensler et al. 2002) and ^~(0.3-0.6)c for Vela (Kargaltsev et 
al. 2003). 

5.3. Comparison of PSR J1119— 6127 to other rotation-powered pulsars and 

magnetars 

There are currently 7 pulsars including PSR J1119— 6127 with magnetar-strength fields 
(B>B c =4.4xl0 13 G, see Table 3), but which have not been identified as magnetars, except 
for PSR J1846— 0258 which most recently revealed itself as a magnetar (Kumar & Safi-Harb 
2008, Gavriil et al. 2008). One would ask then whether all other pulsars belonging to this 
class should be related to magnetars (e.g. by being quiescent magnetars), or whether they 



represent a distinct population of rotation-powered pulsars. To address this open issue, 
we compare the known properties of PSR J1119— 6127 to magnetars and to the 'classical' 
rotation-powered pulsars: 



• Unlike magnetars (even in quiescence) and like the classical rotation-powered pulsars, 
PSR Jl 119— 6127 and its PWN are powered by rotation (L x is much smaller than its 
E) and the pulsar powers a hard PWN. As well its infrared properties resemble more 
the rotation-powered pulsars than the magnetars (Mignani et al. 2007) - the latter are 
known to be more efficient IR emitters. 

• Like magnetars and some rotation-powered pulsars (see more below), J1119— 6127 is de- 
scribed by a two-component (blackbody+power law) model. PSR J1119— 6127's ther- 
mal temperature (A;T~0.21 keV in the BB+PL model and 0.14 keV in the NSA+PL 
model) is intermediate between that of magnetars (£;Tb B ~0.41-0.67 keV, Woods & 
Thompson 2006) and that of the Vela-like, rotation-powered, pulsars with a thermal 
component dominating at energies E<2 keV (see Fig. 5 of Zavlin 2007b). The rela- 
tively high temperature of PSR J1119— 6127 (in comparison with the rotation-powered 
pulsars) can be accounted for by a neutron star cooling model with proton superflu- 
idity in its core and a neutron star mass of ~1.5 M (see Fig. 5 of Zavlin 2007b). Its 
photon index (1.0-2.9 in the BB+PL model and 0.7-2.3 in the NSA+PL model, see 
Table 1) appears relatively hard in comparison to that of most AXPs in quiescence (see 
Table 14.1 of Woods & Thompson 2006) and consistent with that expected from the 
rotation-powered pulsars; but also similar to that of three spectrally-hard magnetars 
in their quiescent state (SGR 1806-20, SGR 1900+14, and AXP IE 1841-045 having 
PL indices of 2.0, 1.0-2.5, and 2.0, respectively; Woods & Thompson 2006). 

• PSR J1119— 6127's X-ray properties are most similar to two young 'Vela-like' pulsars 
with a similar E: PSR B1800-21 and PSR J1357-6429's, rotation-powered pulsars 
with a more 'typical' magnetic field (see Table 3 and references therein). Most no- 
tably, both pulsars power compact and hard PWNe and display an ^ value com- 
parable to PSR J1119— 6127's. Furthermore, like PSR J1119— 6127, the spectrum 
of PSR J 1357— 6429 has a thermal component that dominates the total flux at en- 
ergies below 2 keV (~72% for PSR J1357-642 versus ~87% for PSR J1119-6127) 
and is strongly pulsed despite its lower B-field. Such a high pulsed fraction (>50% 
at E<2 keV, Zavlin 2007b) could be explained by intrinsic anisotropy of the thermal 
emission formed in a magnetized neutron star atmosphere coupled with a strong non- 
uniformity of the surface temperature and magnetic field distributions (V. E. Zavlin 
2008, private communication; Zavlin 2007b). Such an interpretation can also be applied 
to PSR J1119— 6127 whose thermal spectrum is also well fitted with a neutron star 
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atmospheric model (see Table 1 and §5.1; V. E. Zavlin 2008, private communication, 
Zavlin 2007a). 

• Among the high-B pulsars, PSR J1119— 6127 has spin properties (i.e. P, B) most 
similar to those of PSR J1846— 0258 (the Kes 75 pulsar). Both pulsars are also now 
believed to be at about the same distance of ~6 kpc. However PSR J1846— 0258 is 
a non-thermal X-ray pnlsar characterized by a hard power-law ind«@ (1X2) and is 
much more efficient at powering a bright PWN (|(PWN)~2%, see Kumar & Safi-Harb 
2008), so their X-ray properties are very different. Furthermore, PSR J1846— 0258 has 
recently revealed itself as a magnetar as magnetar-like X-ray bursts were detected 
with RXTE and its spectrum softened as its flux brightened by a factor of ~6 (Gavriil 
et al. 2008, Kumar & Safi-Harb 2008). PSR Jl 119-6127 appears to have a sta- 
ble spectrum between 2002 (obsid 2833) and 2004 (obsid 4676) with a count rate of 
~5.5xl0 -3 counts s _1 (§4.1). Obsid 6153 yields a smaller count rate for the pulsar 
(~70% the average count rate) suggesting variability; however if real, we believe that 
it is not significant within error and given that 1) obsid 6153 was relatively short and 
taken the day following obsid 4676, and 2) magnetar outbursts normally cause a much 
brighter enhancement in the X-ray flux. It is possible however that PSR J1119— 6127 
will one day reveal itself as a magnetar after an occasional burst driven by its high 
B-field, just like PSR J1846-0258 recently did. Monitoring PSR Jl 119-6127 and the 
other high-B pulsars in the radio and X-ray bands are needed to address this possibility. 

Given the close resemblance of PSR J1119— 6127 to the other rotation-powered pulsars 
with 'typical' B-fields and a similar E, we conclude that the global X-ray properties of 



PSR J1119— 6127 are not (at least not entirely) determined by the high B value 13 !; and that E 
in particular, and likely the environment which should confine the PWN and therefore would 
affect its morphology and brightness, play a significant role in determining its properties. 
Whether PSR J1119— 6127 is a quiescent magnetar that's currently mainly powered by its 



Finally, our study further suggests that PSR J1119— 6127 is characterized by a highly 
magnetized wind, possibly explaining the evidence for a small torus and prominent jets. The 



12 we note here that the column density towards PSR J1846— 0258 is at least twice as high as that towards 
PSR J1119— 6127, thus possibly hindering the significant detection of any soft blackbody component (Kumar 
& Safi-Harb 2008). 

13 One should also keep in mind that the actual magnetic field strength can be a factor of few off from 
the spin-down value, and so pulsars with comparable B may have substantially different actual fields (see 
Camilo 2008 and references therein.) 




E, but that will one day reveal a magnetar-like identity, remains to be seen. 
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deep search for X-ray emission and PWNs around the other high-B pulsars is warranted to 
address the question whether these pulsars are similarly characterized by highly magnetized 
winds. Given the small E of most known high-B pulsars (Table 3), very long exposures will 
be needed to be sensitive to PWNs of luminosity ^10 _4 -10~ 3 i?. 
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Parameter PL BB+PL NSA+PL a 



N n (10 22 cm" 2 ) 


1 - z -0.3 


i-s±5:i 




1.6±0.2 


r 


3-0±g;i 


^-O.Q 




1.5±0.8 


Norm. b (PL) 


5.7l| 2 x 10- 5 


l-SloJ x 10" 


-5 


8.4+^| x lO" 6 


A;T eff (keV) 




0.21±0.01 




1 4+ - 03 


Norm. b (BB) 




3.5±J;| x 10" 


-6 


1.42xl0 -8 (frozen) 


-^unabs(-P-^) C 


1.7lJ;f x 10~ 13 


6.71^ x 10" 


-14 


6.1±H x 10~ 14 


F unabs c (thermal) 




2-2±g"§ x 10" 


-13 


1.8±i;§ x 10- 13 


L x d 10 33 ergs s" 1 


1 4 +L1 


9 4+0.8 




9 1 +0.8 

^•-"--0.6 


Xl (dof) 


1.356 (31) 


0.948 (29) 




0.977 (30) 



Table 1: Spectral fits to the pulsar. All models are modified by interstellar absorption with 
a column density N H . PL denotes a power-law model, BB denotes a blackbody, and NSA 
refers to the neutron star H atmospheric model in XSPEC version 12.4.0. 



a The NSA model assumes a magnetic field of 10 13 Gauss, a neutron star radius of 10 km and mass of 1.4 
Mq. The distance is also fixed at 8.4 kpc (thus the normalization of 1.4xl0~ 8 ). 

ft PL normalization at 1 keV in units of photons cm~ 2 s _1 kcV -1 . BB normalization in units of L^/D\ a 
where L39 is the luminosity in units of 10 39 ergs s^ 1 and D 10 is the distance in units of 10 kpc. 
c 0. 5-7.0 keV unabsorbed flux (ergs cm~ 2 s _1 ). 
d 0. 5-7.0 keV total luminosity (assuming isotropic emission). 
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Parameter/PWN Region 


Compact PWN 


Southern 'Jet' 


r 


i 1+0.9 
i - ± -0.7 


1 4+0-8 
-■-•^-0.9 


Norm. a 


(2.2+ 2 ;^)xl0* 6 


(3.51^) xlO- 6 


xl (dof) 


1.23 (17) 


1.16 (17) 


Anabs (0.5-7 keV) ergs cm" 2 s" 1 


(1.91^) xlO" 14 


(2.5l?;l)xl0- 14 



Table 2: Power law spectral fits to the PWN. The compact PWN refers to the 6"xl5" 
structure elongated in the north-south direction (see Fig. 4). The southern 'jet' corresponds 
to the new faint and extended jet-like structure (~6"x20") shown in Fig. 5. Nh is frozen to 
1.8 xlO 22 cm -2 . Errors are 2a. 



"Normalization at 1 keV in units of photons cm 2 s 1 kcV 
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Pulsar (PSR) 


E (ergs s v ) 


B (10 13 G) 


PWN? 


L x (psr+pwn) 
F, 


Reference 


J1846-0258 


8.0xl0 36 


4.9 


Y(X) 


0.018-0.0325 a d 2 


Kumar & Safi-Harb (2008) 


J1119-6127 


2.3xl0 36 


4.1 


Y(X) 


5xlO- 4 d 2 4 


this work 


J1734-3333 


5.6xl0 34 


5.2 






Morris et al. (2002) 


J1718-3718 


1.5xl0 33 


7.4 




0.013-4.04 5 


Kaspi k McLaughlin (2005) 


J1814-1744 


4.7xl0 32 


5.5 




<1340 b (i 2 


Pivovaroff et al. (2000) 


J1819-1458 c 


2.9 xlO 32 


5.0 




~12d 2 6 


Reynolds et al. (2006) 


J1847-0130 


1.7xl0 32 


9.4 




<18-47 d d 2 


McLaughlin et al. (2003) 


B1800-21 


2.2xl0 36 


0.4 


Y(X) 


io- 4 d 2 


Kargaltsev et al. (2007) 


J1357-6429 


3.1 xlO 36 


0.8 


Y(X) 


6.5xlO" 5 4 5 


Zavlin et al. (2007b) 



Table 3: The high-magnetic field pulsars sorted in order of decreasing spin-down power. We 
used the ATNF pulsar catalog in compiling the pulsar properties and only selected pulsars 
with S>4xl0 13 G (the QED value) (B is defined in footnote#4, but see footnote#13). The 
last two pulsars, shown only for comparison with PSR J1119— 6127, are rotation-powered 
pulsars with a 'typical' magnetic field, but with E and X-ray properties almost identical to 
PSR Jl 119— 6127. The 4th column refers to the presence (Y) or absence of a PWN in X-rays 
(X). Except for PSR J1846— 0258, none of these pulsars has a radio PWN associated with 
it. The 5th column refers to the ratio of the non-thermal X-ray luminosity of the combined 
pulsar and associated PWN (when detected in X-rays) to the pulsar's spin down power, 
scaled to the distance given in units of d x , where x is the adopted distance to the pulsar in 
units of kpc. 



a the lower (upper) range corresponds to the 2000 (2006) Chandra data; the pulsar has brightened by a factor 
of -6 in 2006. 

fc only an upper limit on the X-ray luminosity was estimated using the spectral properties of IE 2259+586 (an 
AXP) as a template. 

c pulsar classified as a rotating radio transient; the X-ray luminosity of the pulsar is uncertain by an order of 
magnitude. 

d only an upper limit was estimated using AXPs 4U0142+61 or IE 2259+586 as templates. 
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Fig. 1.— 2.0-10.0 keV ASCA image of SNR G292.2-0.5, smoothed with a Gaussian with 
a— 45", with the black circle showing the radio boundary (15' in diameter). The superimposed 
squares mark the locations of the 8'x8' Chandra's S3 chip during AO-3 (right) and AO-5 
(left), covering the western and eastern sides of the remnant, respectively, with the grayscale 
image showing the raw 0.5-10 keV emission seen with Chandra. The central circle marks 
the position of PSR J1119— 6127. North is up and east is to the left. 
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Fig. 2. — The tri-color energy image of PSR J1119— 6127 and associated PWN. The image 
is 2.4'x2.3' in size and is exposure corrected. Red, green, and blue correspond to 0.5-1.15 
keV, 1.15-2.3 keV, and 2.3-10.0 keV, respectively (see §3.1 for details). 



- 23 - 



0.5-3 keV 




S -37'oD' 




ll a 19 m 15" 14 s 13 s 

Right Ascension (J2000) 
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Fig. 3.— The raw soft (0.5-2.0 keV) and hard (2.0-10 keV) images of PSR J1119-6127 and 
associated nebula. These images were not background-subtracted. 



2-10 keV smoothed 




Right Ascension (J2000) 11 19 15 14 13 ' 

Right Ascension (J2000) 

Fig. 4. — The same images as in Fig. 3, adaptively smoothed (with 'csmooth') using a 
Gaussian with a—1" for a significance of detection greater than 5 and up to cr=3" for a 
significance of detection down to 2. The PWN and its elongated morphology in the north- 
south direction is evident in the hard band image. Contours are overlayed on a linear scale, 
chosen to highlight the emission from the jets, and ranging from 2.5 to 17 counts per pixel 
in steps of 2 for the soft band image, and from 2.5 to 7 counts per pixel in steps of 0.75 for 
the hard band image. 
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Fig. 5. — The 2-10 keV image of PSR J1119— 6127 displaying the extended jet-like feature 
southwest of the pulsar. The image has been background-subtracted and smoothed as the 
previous figure. The overlaid contours range from 1.4 to 10 counts per pixel in steps of 1 in 
linear scale and are chosen to highlight the structure of the elongated jet. 
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Fig. 6. — The spectrum of PSR J1119— 6127 fitted with the blackbody model, illustrating 
the excess of emission above ~3 keV unaccounted for by the blackbody model. The lower 
panel shows the ratio of data to plotted model. The 3 colors refer to the 3 observations used 
in our analysis (obsid 2833 in black, 4676 in red, and 6153 in green). 
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Fig. 7. — The spectrum of PSR J1119— 6127 fitted with the BB+PL model shown in Ta- 
ble 1. The lower panel shows the ratio of data to plotted model. The dotted lines show the 
contributions of the blackbody (soft) and power law (hard) components to the spectrum. 
The 3 colors refer to the 3 observations used in our analysis (obsid 2833 in black, 4676 in 
red, and 6153 in green). 



